Goos-Hanchen shifts in frustrated total internal reflection studied with wave packet 

propagation 
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We have investigated that the Goos-Hanchen (GH) shifts in frustrated total internal reflection 
(FTIR) studied with wave packet propagation. In the first-order approximation of the transmission 
coefficient, the GH shift is exactly the expression given by stationary phase method, thus saturates 
an asymptotic constant in two different ways depending on the angle of incidence. Taking account 
into the second-order approximation, the GH shift always depends on the width of the air gap due 
to the modification of the beam width. It is further shown that the GH shift with second-order 
correction increases with decreasing the beam width at the small incidence angles, while for the large 
incidence angles it reveals a strong decrease with decreasing the beam width. These phenomena 
offers the better understanding of the tunneling delay time in FTIR. 

PACS numbers: 42.25.Bs, 42.25.Gy, 73.40.Gk 



It is well known that a light beam totally reflected 
from an interface between two dielectric media under- 
goes lateral shift from the position predicted by geo- 
metrical optics [l|. This phenomenon was referred to 
as Goos-Hanchen (GH) effect and was theoretically ex- 
plained by Artmann's stationary phase method Q and 
Renard's energy flux method Because of the poten- 
tials applications in integrated optics [31, optical waveg- 
uide switch and optical sensors [1,13, the GH shifts 
including other three non-specular effects such as an- 
gular deflection, focal shift and waist-width modifica- 
tion have been extensively investigated in partial reflec- 
tion H i, [H [H El, El], attenuated total reflection 
'14'. T^, and frustrated total internal reflection (FTIR) 
mil?, 18, 19, 20]. 

From a somewhat different perspective, the optical 
tunneling phenomenon in FTIR have attracted much at- 
tention in the last two decades [HI US, [H, 113, HI , be- 
cause of the analogy between FITR and quantum tun- 
neling. Theoretical [H [H and experimental [H [H 
investigations have demonstrated that the GH shifts in 
FTIR play an important role in the supcrluminal tun- 
neling time and the well-known "Hartman effect" [26j . 
which describes that the group delay for quantum parti- 
cles tunneling though a potential barrier saturates to a 
constant for an opaque barrier. Recently, Martinez and 
Polatdemir [23| have studied the effect of width of the 
beam on the GH shift (which is proportional to tunnel- 
ing time) to offer the complementary insights into the 
origin of "Hartman effect" in FTIR. In addition, Haibel 
et al. [1^ once carried out a comprehensive study of 



the GH shift in FTIR as a function of the polarization, 
beam width, and incidence angle in the microwave ex- 
periment, which challenges its theoretical descriptions. 
However, the current expressions of the GH shifts given 
by stationary phase method and energy flux method are 
independent of the beam width. 

The main purpose of this Brief Report is to investi- 
gate that the GH shifts in FTIR by wave packet propa- 
gation. It is shown that the GH shift in the first-order 
approximation of the transmission coefficient is exactly 
the expression of the GH shift given by stationary phase 
method. The GH shift in this case approaches the satu- 
ration value in two different ways depending on the in- 
cidence angle. Taking account into the second-order ap- 
proximation, the GH shift always depends on the width 
of air gap. It is further shown that the GH shift with 
the second-order correction become dependent strongly 
on the beam width. These phenomena offers the better 
understanding of the tunneling delay time in FTIR. 

For simplicity, we consider TE polarized beam inci- 
dent into the double-prism structure with the angular 
frequency to and incidence angle ^o, shown in Fig. [1] 
where a is the width of the air gap. Denote by e, /i and 
n, respectively, the permittivity, permeability and refrac- 
tive index of the prism. For a well-coUimated beam, the 
electric field of the incident beam can be expressed by 
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A(ky)ex:p[i{kxX + kyy)]dky, (1) 
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where — nk cos 9, ky — nk sin 9, k = lo/c, n ^ \/e/i, c 
is the speed of light in vacuum, is the incident angle of 
the plane wave component under consideration, and time 
dependence exp(— zcji) is implied and suppressed. For a 
Gaussian-shaped incident beam whose peak is assumed 
to be located at a; = 0, 
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FIG. 1: Schematic diagram of the GH shift in FTIR, where 
the width of the air gap is a, St represents the GH shift of the 
transmitted beam. 



= 0,y) = exp ^-^j exp(i/cyoy), (2) 

its angular spectral distribution is also a Gaussian func- 
tion, A{ky) = Wyex'p[—{wy/2){ky — kyo)'^], around its 
central kyo — k sin^Oi Wy — wq/ cos^O) Wq is the width of 
the beam at waist. According to Maxwell equations and 
boundary conditions, the field of the transmitted beam 
is found to be 
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TA{ky) exp{i[fc-c(.T — a) + kyy]}dky, 

(3) 

with the transmitted coefficient T = eKp(i(f))/f is given 
by 

/exp(i(/)) — coshKa-|- — [ — — — J sinh/ta, (4) 
2 y /v kx J 

where k — {ky ~ k^Y^^. 

Firstly, we look at the GH shift in the first-order ap- 
proximation of the transmission coefficient. Expand the 
exponent in Taylor series at kyo, and retain up to the 
first-order term, then we will obtain 



T[ky) — exp[lnT(A;y)] « Tq exp 



1 dT 

Tq dkyo 



(ky kyo) 



(5) 

where Tq = T{kyo), and d/dkyo denotes the derivative 
with respect to ky evaluated at ky = kyQ. Introduce two 
real parameters and L" defined as. 



L = L't+ iL'l = 



dT 



Tq dk 



yO 



(6) 



then, in terms of the phase and magnitude of T, we will 
have 



dk 



yO 



and 



dk, 



yO 



lli\T{ky)\. 



Substituting Eq. ([5]) into Eq. ([3]) and employing the 
paraxial approximation condition, 

(k — k o)^ 

kx ~ kxQ - {ky - kyo) tan^o - — | (7) 

ZiK, COS i/Q 

we obtain the transmitted beam at x = a, 
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(8) 

It is clear that the lateral shift L'^ — —d(j)/dkyQ is the same 
as the one obtained by the stationary phase method [3], 
and is given by 
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(9) 
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'^xQ 



sinh 2^00 
2KQa 



where Sc = akyQ/nQ. When the width of the air gap is 
large enough, that is, a 3> the GH shift will tends 
to a constant. 



asymp ' 



lim Sj 

KQa — *oo 
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yO 



KokxO 



(10) 



With increasing the air gap the GH shift reaches a asymp- 
totic constant, which is in agreement with the experimen- 
tal results , and is also closely related to the counter- 
intuitive "Hartman effect" of the tunneling delay time in 
the limit of an opaque barrier [H, [23| . 

More interestingly, what we emphasized here is that 
the GH shift approaches the saturation value in two dif- 
ferent ways depending on the angle of incidence. The GH 
shift ([9]) can be expressed by the following form: 
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(11) 



where go = {k^Q — Kq) tanh KQa/2kxQKQ. Keeping the next 
to leading term for large a shows the approach to asymp- 
totic value by 



asymp 
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(kq -I- k^f^y 



-2KQa 



(12) 



It is clearly evident from the above expression that for 
Kg — fc^Q < the GH shift increases monotonically to 
reach the saturation value, while for Kq — k^Q > it 
reaches the saturation value from above, that is, there is 
a hump before it attains saturation. Therefore, when the 
necessary condition for incident angle. 



iQ > Op = sm 



1+7^2 



(13) 



'P-— V 2^2 ' 
is satisfied, the GH shift can approach the saturation 
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FIG. 2; (Color online) Dependence of the GH shifts (in unit 
of A) on the air gap width a in FTIR, where A — 32.8mm and 
n — 1.605, the solid, dotted and dashed curves correspond to 
the incidence angles 80°, 75°, and 70°, respectively. 
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FIG. 3: (Color online) Dependence of GH shift (in unit of 
A) on width of air gap with different beam widths, where (a) 
9o = 45° (b) 6o — 75° , and other parameters are the same as 
in Fig. [5] The solid corresponds to GH shift in the first-order 
approximation, the dashed and dotted curves correspond to 
the GH shifts in the second-order approximation. 



limit with negative slope and can be larger than the sat- 
uration value for the intermediate values of the air gap. 
The position of the maximum of the hump (ao) is given 
by 



then, with the phase and magnitude of T{ky). we have 
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ao 



3/c^o' 



2ko (fc^Q + kI) 



2\2 



(14) 



and 



a/tj^O 



InlTl 



Substituting expression into Eq. ^ using paraxial 
condition ([7]), and neglecting some unimportant factors, 
we finally obtain the following field of the transmitted 
beam at a; = a. 



*t(a,?/) « Toexp 



exp 



1 



Fig. [5] shows that for the large a the GH shift is 
independent of the width a of the air gap hence it sat- 
urates a asymptotic constant, where A — 32.8mm and 
n = 1.605 (corresponding to critical angle 9c = 38.5° for 
total reflection and 9p = 56.4°) [l^. Furthermore, the 
GH shift approaches the asymptotic limit from above for 
6o > 9p and from below for 6q < Op. These phenomenon 
is not due to the interference time |27|, and does result 
from the interference between the incident and reflected 
beams. Of course, it can also been seen from the rela- 
tionship between the GH shift and group delay discussed 
in Ref. [2^ that the delay time in FTIR also saturates to 

a constant from above for 6q > 0p 29] in the same way where r]t — L'/ /wty, wtf — (z 



as that in the quantum tunneling for < Vq /2 [30|, |31| 
since the self-interference delay time that comes from the 
overlap of incident and reflected waves in front of barrier 
is of great importance j33 | . 

Next, what as follows we will show the influence of the 
beam waist width on the GH shift in FTIR. To this end, 
we consider the exponent of the transmission coefficient 
is approximated to the second-order term. 



y- L't 



run 
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(17) 
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and 



"ty 



Wy + Fl' correspond to the angular deflection, focal shift 
and waist-width modification, respectively [Ij]. The GH 
shift in this case can be expressed by, 
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T{ky) w Toexp 
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2 dkyQ \ T dk 



1 dT 



To dkyQ 



(15) 



Introducing two new real parameters F^ and F^ defined 

as 



F^F + iF, 



1 dT 



dkyQ \ T dky 



(16) 



Obviously, the second term on the right-handed side of 
Eq. (fT8|) is a second-order correction, which leads to 
the dependence of the GH shift on the beam width. In 
addition, it also results in the dependence of the GH shift 
on the width of air gap in the opaque barrier limit. 

Fig. [3] demonstrates that the GH shift in the second- 
order approximation depends on the beam width, where 
(a) 9q — 45° (b) 9q = 75°, and other parameters are 
the same as in Fig. [51 Compared with Fig. [2] discussed 
above, the GH shift becomes dependent on the width a 
in the limit of an opaque barrier, due to the second cor- 
rection. When the beam width is large, that is, the di- 
vergence angle becomes small, the correction to GH shift 
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can be neglected, thus for a well-coUimated beam the 
GH shift is in agreement with that given by the station- 
ary phase method. More importantly, Fig. [3] shows that 
the GH shift increases with decreasing the beam width at 
6*0 = 45°, while the GH shift for 6*0 = 75° shows a strong 
decrease with decreasing the beam width. As shown in 
Fig. [3] (a) , the GH shift becomes dependent linearly on 
the width of air gap, because the Fourier components of 
the incident beam above the critical angle are strongly 
depressed so that the plane wave components just be- 
low the critical angle start to dominate. That is to say, 
when the incidence angle is larger than but close to the 
critical angle, the wave vector filter is more pronounced 
for a larger beam width, the transmission is essentially 
not tunneling at all, thus the GH shift increases with in- 
creasing the width a, as one expects classically. This also 
implies the violation of "Hartman effect" for the quan- 
tum tunneling in time domain [s^l • 

Finally, we have a brief look at the microwave exper- 
iment on GH shifts. It was once argued [l^ that the 
influences of beam width and incidence angle challenge 
the current descriptions of the GH shift in FTIR. Fig. 
13] (a) shows the GH shift increases with deceasing beam 
dimension corresponding to the beam waist width, which 
is agreement with the experimental results in Ref. [l^ 
where the physical parameters are the same as those in 
Fig. [51 In addition, it is also predicted in Fig. [3] (b) that 
the GH shift will decreases with deceasing beam waist, 
when the incidence angle is away from the critical an- 



gle. In a word, the improved formula of GH shift (fTS)) 
with the modification of beam width can give better un- 
derstanding of the GH shift in FTIR theoretically and 
experimentally. 

To summary, we have investigated the GH shifts in 
FTIR by wave packet propagation. It is found that the 
GH shift in the first-order approximation of the transmis- 
sion coefficient, which is exactly the expression of the GH 
shift obtained by stationary phase method, approaches 
the saturation value in two different ways depending on 
the angle of incidence. The explicit expression of the GH 
shift in the second-order approximation shows the strong 
dependence on the beam width. It is further shown that 
the GH shift with the second-order correction increases 
with decreasing the beam width at the small incident 
angles, while for the large incident angles the GH shift 
reveals a decrease with decreasing the beam width. All 
these theoretical results can be applicable to explain the 
experiment on GH shifts [19] and offer a hint to the better 
understanding of tunneling delay time in FTIR [2^ . 
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